home *** CD-ROM | disk | FTP | other *** search
/ Packard Bell - Internet on a CD / internet on a cd.cdr / Internet / sites / Clementine_NASA / clemdcmp.c < prev    next >
Encoding:
C/C++ Source or Header  |  1998-07-16  |  17.2 KB  |  528 lines
  1. /*************************************************************************     
  2. *                   ----------------
  3. *                   |  CLEMDCMP    |
  4. *                   ----------------
  5. *
  6. *_TITLE CLEMDCMP Decompress clementine image into various formats
  7. * 
  8. *_DESC  CLEMDCMP will decompress a clementine image and format it into 
  9. *       one of three possible formats:
  10. *       1) decompressed pds image file, with pds labels, the historgram
  11. *           object, and an image object, either the browse image or
  12. *           the full image
  13. *       2) decompressed image file, no labels
  14. *       3) decompressed gif image
  15. *       4) decompressed tiff image
  16. *
  17. *       CLEMDCMP has four command-line arguments.  The first argument
  18. *       indicates the output format of the decompressed image and
  19. *       the following are the options:
  20. *         -p = decompressed pds image file
  21. *         -n = decompressed unlabeled file
  22. *         -t = decompressed tiff file
  23. *         -g = decompressed gif file
  24. *       The second argument indicates which image object to output:
  25. *         -i = full resolution image
  26. *         -b = browse image
  27. *       The third argument is the filename of the input compressed clementine 
  28. *         data file.
  29. *       The fourth argument is the filename of the output decompressed
  30. *         clementine image.
  31. *
  32. *_HIST  Aug 27 1994  Tracie Sucharski, USGS, Flagstaff Original Version
  33. *
  34. **************************************************************************/
  35.  
  36. #include <stdio.h>
  37. #include <stdlib.h>
  38. #include <string.h>
  39.  
  40. #include "pds.h"
  41.  
  42. #define MM_TYPE 0x4D4D
  43. #define II_TYPE 0x4949
  44. #define BLKSIZE 32768L
  45.  
  46. main(int argc, char *argv[])
  47. {
  48.   int c;                         /* character for command-line arguments */
  49.   char format;                   /* Output file format */
  50.   char img;                      /* Output image or browse image */
  51.   char infile[128];              /* Input file name */
  52.   char outfile[128];             /* Output file name */
  53.   long lines;                    /* Number of lines and sample of image */
  54.   long samps;
  55.   int  i;                       
  56.   FILE *fpi, *fpo;               /* Input and output file pointers */
  57.   int ret;                       /* Return code */
  58.   char b[2];
  59.   short int *j;
  60.   char bitord;
  61.   unsigned int k, nblocks, rem;
  62.   
  63.   PDSINFO *p;
  64.  
  65.   int labmod(char *text, char img, char bitord);
  66.   int writetif(FILE *fp, long nl, long ns, CHARH *buf, char bitord);
  67.  
  68.  
  69.  
  70.   for (i = 1; i < argc; i++) {
  71.     if (argv[i][0] == '-') {
  72.     argv[i][1] = tolower(argv[i][1]);
  73.     switch (argv[i][1]) {
  74.       case 'p': case 'n': case 'g': case 't':
  75.     format = argv[i][1];
  76.     break;
  77.       case 'i': case 'b':
  78.     img = argv[i][1];
  79.     break;
  80.       default:
  81.     printf("ERROR:CLEMDCMP-Illegal option chosen %c\n", c);
  82.     exit(1);
  83.       }
  84.   }
  85.     else
  86.       break;
  87.   }
  88.  
  89.   if (argc-i != 2) {
  90.     printf("ERROR:CLEMDCMP-You must specify an input and output file\n");
  91.     exit(1);
  92.   }
  93.  
  94. /*************************************************************************
  95. *  Default is to output the full image in pds format
  96. *************************************************************************/
  97.   if (format == '\0') format = 'p';
  98.   if (img == '\0') img = 'i';
  99.  
  100.  
  101.   strcpy(infile, argv[i]);
  102.   strcpy(outfile, argv[i+1]);
  103.  
  104.  
  105. /************************************************************************
  106. *  Determine the bit order
  107. ************************************************************************/
  108.  
  109.   j = (short int *) b;
  110.   b[0] = 1;
  111.   b[1] = 0;
  112.   if (*j == 1) bitord = 'l';
  113.   else bitord = 'm';
  114.  
  115.  
  116. /****************************************************************************
  117. *  The following call will open the compressed Clementine image and create
  118. *  a structure containing the file objects.
  119. ****************************************************************************/
  120.   p = PDSR(infile, &lines, &samps);
  121.  
  122. /****************************************************************************
  123. *  If user chooses to output the browse image, make sure it exists.
  124. ****************************************************************************/
  125.   if (img == 'b' && !(p->brw_imag)) {
  126.     printf("ERROR:CLEMDCMP-This file does not contain a browse image");
  127.     exit(1);
  128.   }
  129.  
  130. /******************************************************************************
  131. *  Open output file
  132. ******************************************************************************/
  133.   fpo = fopen(outfile, "wb");
  134.  
  135. /*****************************************************************************
  136. *  For the pds formatted option, the labels will need to be modified before 
  137. *  writing out.  The uncompressed output file will either be the image or 
  138. *  the browse image, so the pointer to the browse image will be deleted, 
  139. *  and the image pointer will be modified to point to the start of the data.  
  140. *  The pointer to the histogram will also need to be modified.  The 
  141. *  ENCODING_COMPRESSION_RATIO will no longer be applicable, so the value is 
  142. *  changed to "N/A". 
  143. ******************************************************************************/
  144.   if (format == 'p') {
  145.     printf("Writing PDS labeled file.\n");
  146.     if (labmod(p->text, img, bitord) < 0) {
  147.       printf("ERROR:CLEMDCMP-Error reformatting labels");
  148.       exit(1);
  149.     }
  150.     fwrite(p->text, strlen(p->text), 1, fpo);
  151.   }
  152.  
  153. /*****************************************************************************
  154. *  If the output is a pds file write the histogram data followed by the image 
  155. *  data, either the browse image or the full image.  If the output is an
  156. *  unlabeled image file,  write only the image data out.
  157. *****************************************************************************/
  158.   if (format == 'p') fwrite(p->hist, sizeof(long), 256, fpo);
  159.  
  160.   if (format=='p' || format=='n') {
  161.     if (format=='n') printf("Writing raw unlabeled file.\n");
  162.  
  163.     if (img == 'i') {
  164. #ifdef __TURBOC__
  165.       nblocks = (lines*samps) / BLKSIZE;
  166.       rem = (lines*samps) % BLKSIZE;
  167.       for (k=0; k<nblocks; k++, p->image+=BLKSIZE)
  168.     fwrite(p->image,1,BLKSIZE,fpo);
  169.       if ( rem > 0 )
  170.     fwrite(p->image,1,rem,fpo);
  171. #else
  172.       fwrite(p->image, lines*samps, 1, fpo);
  173. #endif
  174.       if (format=='p')printf("PDS labeled file was successfully written.\n");
  175.       if (format=='n') {
  176.         printf("Raw unlabeled file was successfully written.\n");
  177.     printf("\nThe size of the raw image created is %ld lines "
  178.                "by %ld samples with no\n header or label data.\n\n", 
  179.                lines, samps);
  180.       }
  181.     }
  182.  
  183.     if (img == 'b') {
  184.       fwrite(p->brw_imag, p->browse_nrows*p->browse_ncols, 1, fpo);
  185.  
  186.       if (format=='p')printf("PDS labeled file was successfully written.\n");
  187.       if (format=='n') {
  188.     printf("Raw unlabeled file was successfully written..\n");
  189.     printf("\n\nThe size of the raw image created is %ld lines "
  190.            "by %ld samples with no\n header or label data.\n\n", 
  191.            p->browse_nrows, p->browse_ncols);
  192.       }
  193.     }
  194.  
  195.   }
  196.  
  197.  
  198.  
  199.   if (format == 't') {    /*  TIFF format  */
  200.     printf("Writing TIF file.\n");
  201.     if (img == 'i') {
  202.       if (writetif(fpo, lines, samps, p->image, bitord) < 0) {
  203.     printf("ERROR:CLEMDCMP-Error writing full image in tiff format");
  204.     exit(1);
  205.       }
  206.     }
  207.  
  208.     if (img == 'b') {  /* Write browse image */
  209.       if (writetif(fpo, p->browse_nrows, p->browse_ncols, p->brw_imag, bitord) < 0) {
  210.     printf("ERROR:CLEMDCMP-Error writing browse image in tiff format");
  211.     exit(1);
  212.       }
  213.     }
  214.     printf("TIF file was successfully written.\n");
  215.   }
  216.  
  217.   if (format == 'g')  {   /*  GIF format  */
  218.     if (img == 'i') {
  219.       if (writegif(fpo, lines, samps, p->image) < 0) {
  220.     printf("ERROR:CLEMDCMP-Error writing full image in gif format");
  221.     exit(1);
  222.       }
  223.     }
  224.  
  225.     if (img == 'b') {
  226.       if (writegif(fpo, p->browse_nrows, p->browse_ncols, p->brw_imag) < 0) {
  227.     printf("ERROR:CLEMDCMP-Error writing browse image in gif format");
  228.     exit(1);
  229.       }
  230.     }
  231.   }    
  232.  
  233.   
  234.   fclose(fpo);
  235.   exit(0);
  236. }
  237.  
  238.  
  239.  
  240.  
  241.  
  242. int labmod(char *text, char img, char bitord)
  243. /*****************************************************************************
  244. *
  245. *_TITLE  LABMOD - modify the clementine label to reflect the decompression
  246. *
  247. *_ARGS    Type    Variable     I/O  Description
  248. *_PARM    char    *text;        I   Pointer to clementine label
  249. *_PARM    char    img           I   Image Type (Full image or browse image)
  250. *_PARM    char    bitord        I   Bit order of current machine
  251. *_PARM    int     *ret          O   Return code
  252. *                                   0 - OK
  253. *
  254. *_DESC  LABMOD will modify the clementine label to reflect the decompression,
  255. *       and the possible rearranging of the browse image and image object.
  256. *
  257. *_HIST  Apr 13 1994 Tracie Sucharski, USGS, Flagstaff Original Version
  258. *       Jun 27 1994 Tracie Sucharski,  Fixed bug when writing out the
  259. *                       browse image.
  260. *       Jul 14 1994 Tracie Sucharski,  Changed output value of ENCODING_
  261. *                       TYPE keyword.
  262. *       Aug 23 1994 Tracie Sucharski,  Added NOTE to image object if the
  263. *                       output image is the browse image,  also correct
  264. *                       the DATA_TYPE keyword in the histogram object
  265. *                       to indicate the correct byte order.
  266. *
  267. ****************************************************************************/
  268. {
  269.     int lbllen;                     /* Length of incoming label */
  270.     int nlbllen;                    /* Length of outgoing label */
  271.     char *start, *end;              /* Index pointers */
  272.     char sdummy[17];
  273.     int nc;                         /* Number of characters */
  274.     char byte[4];                   /* Starting bytes of objects */
  275.     int hbyte;                      /* Starting byte of Image Historgram */
  276.     int ibyte;                      /* Starting byte of Image */
  277.     int bbyte;                      /* Starting byte of Browse Image */
  278.     int diff;                       /* Difference in bytes between input */
  279.                                     /* label and output label */
  280.     char *labels;                   /* Temporary buffer to hold labels */
  281.  
  282.  
  283.     lbllen = strlen(text);
  284.     labels = (char *)malloc(lbllen);
  285.  
  286.     start = strstr(text, "^IMAGE_HISTOGRAM ");
  287.     sscanf(start, "%s = %d", sdummy, &hbyte);
  288.  
  289.     start = strstr(text, "^IMAGE ");
  290.     sscanf(start, "%s = %d", sdummy, &ibyte);
  291.  
  292.     start = strstr(text, "^BROWSE_IMAGE ");
  293.     sscanf(start, "%s = %d", sdummy, &bbyte);
  294.  
  295. /****************************************************************************
  296. *  Get rid of pointer to browse image.
  297. ****************************************************************************/
  298.     end = strchr(start,'\n');
  299.     strcpy(start, end+1);
  300.  
  301. /***************************************************************************
  302. *  If browse image is being written, get rid of IMAGE object, rename
  303. *  BROWSE_IMAGE object to IMAGE and add a note to indicate the output
  304. *  image is the browse image.
  305. ***************************************************************************/
  306.     if (img == 'b') {
  307.       start = strstr(text, "OBJECT = IMAGE\n");
  308.       end = strstr(start, "END_OBJECT");
  309.       strcpy(start, end+11);
  310.  
  311.       start = strstr(text, " BROWSE_IMAGE");
  312.       end = strchr(start, '\n');
  313.       strncpy(labels, text, (start-text)+1);
  314.       *(labels+(start-text+1)) = '\0';
  315.       strcat(labels, "IMAGE");
  316.       strcat(labels, end);
  317.       strcpy(text, labels);
  318.  
  319.       end = strstr(start, "END_OBJECT");
  320.       strncpy(labels, text, (end-text)+1);
  321.       *(labels+(end-text)) = '\0';
  322.       strcat(labels, "  NOTE          = \"Averaged subsampled EDR image\"\n");
  323.       strcat(labels, end);
  324.       strcpy(text, labels);
  325.  
  326.     }
  327.  
  328.     if (img == 'i') {
  329.  
  330. /**************************************************************************
  331. *  Get rid of the BROWSE_IMAGE object.
  332. ***************************************************************************/
  333.       start = strstr(text, "OBJECT = BROWSE_IMAGE");
  334.       end = strstr(start, "END_OBJECT");
  335.       strcpy(start, end+11);
  336.  
  337. /***************************************************************************
  338. *  If writing IMAGE object change ENCODING_COMPRESSION_RATIO and
  339. *  ENCODING_TYPE keyvalues to "N/A", since the output file is uncompressed.
  340. ***************************************************************************/
  341.       start = strstr(text, "ENCODING_TYPE ");
  342.       end = strchr(start, '\n');
  343.       strncpy(labels, text, (end-text)-1);
  344.       *(labels+(end-text-1)) = '\0';
  345.       strcat(labels, " DECOMPRESSED\"");
  346.       strcat(labels, end);
  347.       strcpy(text, labels);
  348.  
  349.       start = strstr(text, "ENCODING_COMPRESSION_RATIO ");
  350.       end = strchr(start, '\n');
  351.       strncpy(labels, text, (start-text)+29);
  352.       *(labels+(start-text+29)) = '\0';
  353.       strcat(labels, "\"N/A\"");
  354.       strcat(labels, end);
  355.       strcpy(text, labels);
  356.  
  357.     }
  358.  
  359. /****************************************************************************
  360. *  If the bitorder is MSB which is different from that in the label(LSB),
  361. *  change the label.
  362. ****************************************************************************/
  363.     if (bitord == 'm') {
  364.       start = strstr(text, "LSB_INTEGER");
  365.       *start = 'M';
  366.     }
  367.  
  368. /****************************************************************************
  369. *  Adjust the object pointers.
  370. ****************************************************************************/
  371.     nlbllen = strlen(text);
  372.     hbyte = nlbllen + 1;
  373.     ibyte = hbyte + 1024;
  374.  
  375. /*****************************************************************************
  376. *  Now that there are new pointer values, write them out to the labels.
  377. *****************************************************************************/
  378.     sprintf(byte, "%d", hbyte);
  379.     nc = strlen(byte);
  380.     start = strstr(text, "^IMAGE_HISTOGRAM ");
  381.     strncpy(start+19, byte, nc);
  382.  
  383.     sprintf(byte, "%d", ibyte);
  384.     nc = strlen(byte);
  385.     start = strstr(text, "^IMAGE ");
  386.     strncpy(start+19, byte, nc);
  387.  
  388.     free(labels);
  389.     return(0);
  390.  
  391. }
  392.  
  393.     
  394.  
  395. /*************************************************************************
  396. *                       ------------
  397. *                       | WRITETIF |
  398. *                       ------------
  399. *
  400. *_TITLE  WRITETIF  Writes decompressed Clementine data into TIF format
  401. *
  402. *_DESC   WRITETIF takes decompressed Clementine data and writes an 
  403. *        uncompressed TIF formatted image.
  404. *
  405. *_HIST   May 04 1994  Tracie Sucharski, USGS, Flagstaff Original Version
  406. *
  407. *************************************************************************/
  408.  
  409. int writetif(FILE *fp, long nl, long ns, CHARH *buf, char bitord)
  410. {
  411.   unsigned int j, nblocks, rem; /* Added by Luis Perez 06/29/94 */
  412.  
  413.   int fputword(FILE *fp, short int n);
  414.   int fputlong(FILE *fp, long n);
  415.  
  416.  
  417. /********************************************************************
  418. *  Write out the TIF header
  419. ********************************************************************/
  420.   if (bitord == 'm')  fputword(fp, MM_TYPE);
  421.   if (bitord == 'l')  fputword(fp, II_TYPE);
  422.   fputword(fp,42);
  423.   fputlong(fp,8L);
  424.  
  425. /********************************************************************
  426. *  Construct the Image File Directory (IFD)
  427. ********************************************************************/
  428.   fputword(fp, 8);     /* Eight Tags */
  429.  
  430.   fputword(fp, 254);   /* NewSubfileType */
  431.   fputword(fp, 4);     /* Long */
  432.   fputlong(fp, 1L);
  433.   fputlong(fp, 0L);
  434.  
  435.   fputword(fp, 256);   /* ImageWidth */
  436.   fputword(fp, 3);     /* Short */
  437.   fputlong(fp, 1L);    
  438.   if (bitord == 'm')
  439.     fputlong(fp, ns*65536); /* Shift value into the two high order bytes */
  440.   if (bitord == 'l')
  441.     fputlong(fp, ns);
  442.  
  443.   fputword(fp, 257);   /* ImageLength */
  444.   fputword(fp, 3);     /* Short */
  445.   fputlong(fp, 1L);
  446.   if (bitord == 'm')
  447.     fputlong(fp, nl*65536);  /* Shift value into the two high order bytes */
  448.   if (bitord == 'l')
  449.     fputlong(fp, nl);
  450.  
  451.   fputword(fp, 258);   /* BitsPerSample */
  452.   fputword(fp, 3);     /* Short */
  453.   fputlong(fp, 1L);
  454.   if (bitord == 'm') 
  455.     fputlong(fp, 524288L); /* 8 Shifted into the two high order bytes */
  456.   if (bitord == 'l')
  457.     fputlong(fp, 8L);
  458.  
  459.   fputword(fp, 259);   /* Compression */
  460.   fputword(fp, 3);     /* Short */
  461.   fputlong(fp, 1L);
  462.   if (bitord == 'm')
  463.     fputlong(fp, 65536L);    /* No compression */
  464.   if (bitord == 'l')
  465.     fputlong(fp, 1L);
  466.  
  467.   fputword(fp, 262);   /* PhotometricInterpretation */
  468.   fputword(fp, 3);     /* Short */
  469.   fputlong(fp, 1L);
  470.   if (bitord == 'm')
  471.     fputlong(fp, 65536L);  /* 1 Shifted into the two high order bytes */
  472.   if (bitord == 'l')
  473.     fputlong(fp, 1L);
  474.  
  475.   fputword(fp, 273);   /* StripOffsets - Start of image data */
  476.   fputword(fp, 4);     /* Long */
  477.   fputlong(fp, 1L);
  478.   fputlong(fp, 110L);
  479.  
  480.   fputword(fp, 277);   /* SamplesPerPixel */
  481.   fputword(fp, 3);     /* Short */
  482.   fputlong(fp, 1L);
  483.   if (bitord == 'm')
  484.     fputlong(fp, 65536L);
  485.   if (bitord == 'l')
  486.     fputlong(fp, 1L);
  487.  
  488.   fputlong(fp, 0L);
  489.  
  490. /***********************************************************************
  491. *  Write image data
  492. ***********************************************************************/
  493.   /*!!!!!! BE CAREFUL HERE !!!!!!*/
  494.   /* The third argument in fwrite is of type size_t, which under Turbo C
  495.      becomes an unsigned short int (16 bit integer). This means you have
  496.      to write the image data out in blocks of 64K bytes or less (32K is a
  497.      good number).
  498.      Added by Luis Perez 06/29/94. */
  499. #ifdef __TURBOC__
  500.   nblocks = (nl*ns) / BLKSIZE;
  501.   rem = (nl*ns) % BLKSIZE;
  502.   for (j=0; j<nblocks; j++, buf+=BLKSIZE)
  503.      fwrite(buf,1,BLKSIZE,fp);
  504.   if ( rem > 0 )
  505.      fwrite(buf,1,rem,fp);
  506. #else
  507.   fwrite(buf, 1, nl*ns, fp);
  508. #endif
  509.  
  510.   return(0);
  511.  
  512. }
  513.  
  514.  
  515.  
  516.  
  517. int fputword(FILE *fp, short int n)
  518. {
  519.   fwrite(&n, 2, 1,fp);
  520.   
  521. }
  522.  
  523. int fputlong(FILE *fp, long n)
  524. {
  525.   fwrite(&n, 4, 1, fp);
  526.  
  527. }
  528.